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Abstract 

We study the acoustic scattering properties of a phononic crystal designed to behave as a gradient 
index lens in water, both experimentally and theoretically. The gradient index lens is designed 
using a square lattice of stainless-steel cylinders based on a multiple scattering approach in the 
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homogenization limit. We experimentally demonstrate that the lens follows the graded index 
equations derived for optics by mapping the pressure intensity generated from a spherical source at 
20 kHz. We find good agreement between the experimental result and theoretical modeling based 
i on multiple scattering theory. 
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Composed of ordered arrays of scatterers similar to atoms in a conventional solid, 
phononic crystals (PnC) are a class of metamaterial designed to control acoustic wave 
propagation in a medium. PnCs have been proposed for a broad range of applications in 
wave acoustics, with acoustic lensing^^ featuring prominently in the literature due in part 
to the ease with which focusing can be achieved by altering a crystal's gemoetric shape^ 
or compositional structure.^® Although negative index lenses have received much attention 
due to their potential for near-field imaging,^! some positive index solutions such as the 
acoustic analogue of the optical graded index lens^S^ have not yet been explored experi- 
mentally. In addition, the majority of PnC experiments have been performed in air] 1 * 6 * 11 * ^ 
where the large density contrast with respect to the constituent scatterers in the crystal 
(typically metals) allows the scatterers to be treated as rigid. We demonstrate below that 
despite the physical limitation in impedance contrasts between an aqueous medium and the 
scattering elements, it is possible to design a PnC that behaves as an ideal graded index 
lens (GIL) in water based on a fully elastic multiple scattering theory (MST)P 12 * 13 ^ 

Figure []Ja) shows a plan view schematic of the GIL design. The axes of Fig. [l](a) and 
throughout the paper are oriented with the lens center at position (x,y) = (0,0). The GIL 
is made up of 75 stainless steel cylinders (T-316) that are 75 cm in length and arranged in 
a square lattice with spacing a = 1.8 cm and dimensions 5a x 15a. Figure [Tfb) plots the 
cylinder radii R(y), which are stepped toward zero at each successive layer above and below 
the central axis (y = 0) of the GIL. In the homogenization limit (propagation wavelength 
A > 4a) each stratified layer in Fig. [T^a) can be treated as an effective medium. MST^ 
is used to calculate each layer's effective sound speed c e //, which is inversely proportional 
to the filling fraction of the cylinders. The layers will have an effective refractive index 
n e// — Cb/ceff (o> — 1470 hi/s is the sound speed in water) that is maximal at the center of 
the GIL and decreases to that of water at the edges. Our choice of R(y) in Fig. [jjb) produces 
a graded n e ff that obeys the same relation as an optical GILp^l n e // = n (l — a 2 ?/ 2 ) 1 / 2 , 
where hq is the refractive index at the central layer. Our design results in no = 1.2 and 
a = 0.04 cm" 1 . 

Figure [l](c) shows an image of the GIL submerged in a 6 x 6 x 4 m 3 isolation tank. The 
cylinders are mounted between reinforced Plexiglas plates to provide stability. Acoustic 
waves are produced by a 10 cm-diameter spherical source at 20 kHz (A ~ 4a). The wave 
propagation is measured in the time-domain using monitoring hydrophones at a sampling 
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rate of 1 MHz. Hydrophones are mounted to the source and onto a translational 3-axis 
Velmex VXM® positioning system. The transmission intensity is measured by averaging 
over a 10-cycle pulse from the source; this pulse is long enough to approximate a continuous 
wave measurement, while being short enough to prevent contamination from reflections off 
the surfaces of the tank. We have experimentally verified that the intensity Pq produced by 
the source in the absence of the GIL drops radially in proportion to 1/r 2 . 

Figure ^a) shows the normalized pressure amplitude P/Pq measured after transmission 
through the GIL on the side opposite the source (x > 0). The source is located at (x,y) = 
(—196,0) cm, and both the source and the translational hydrophone are positioned in the 
plane bisecting the axial center of the cylinders. The GIL is shown schematically to scale and 
at its proper location in each figure throughout the paper. The data in Fig. |2](a) is measured 
2.144 ms after the initial cycle began to leave the source. This time gives a snapshot when 
the pulse is centered on an enhancement in signal amplitude observed in the vicinity of 
x ~ 80 cm. 

Figure [2](b) shows the normalized intensity averaged over the 10-cycle pulse and obtained 
from the same data set shown in Fig. |2]^a). As in Fig. |2](a), a clear focusing peak is observed 
centered close to x ~ 80 cm. In Fig. [2^c) we show a two-dimensional MST calculation^^ of 
the total pressure intensity (incident plus scattered) derived by placing a continuous-wave 
cylindrical source at (x,y) = (—196,0) cm. The calculation assumes the cylinders to be 
a penetrable elastic.^^ As with the experimental data, the simulated pressure intensity is 
normalized to that of the source in the absence of the GIL. The source amplitude is a Hankel 
function P = H^\kr) with wavevector k = n /2a. The MST simulation also shows a clear 
focusing peak, but with two important differences: (1) the measured intensity is ~ 2 times 
larger than the simulation, and (2) the simulated focusing peak is slightly farther from the 
GIL and decays more slowly. 

To quantify whether our GIL design behaves as an ideal lens, in Fig. [3]^a) we present 
measurements of the focusing peak along the central axis of the lens (y = 0) for different 
source positions d s . For each d s , a large- amplitude peak is observed above x > 60 cm, while 
smaller peaks are also observed closer to the GIL. As the source is moved closer to the GIL, 
the large-amplitude peak moves away in qualitative agreement with the expected behavior 
of a lens. In Fig. |3](b) we show MST calculations along y = for source positions similar 
to those in the experiment. On initial inspection it appears that the theory shows slowly 
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decaying focusing peaks that change very little with d s . However, expansion of the region 
around the focusing peaks [Fig. ^b) inset] reveals that the peak positions move away in a 
manner similar to the experiment. 

We now analyze the experimental data in Fig. |3^a) above x > 62 cm to determine whether 
the focusing positions in this region follow the ideal lens equations. For an ideal lens, the 
focusing peak positions d p should scale with d s as l/d p = 1/f — l/d s . The focal length / of 
a GIL can be approximated asp^l 

/« (!) 

HqCy sin at 

where t = 5a is the thickness of the GIL. Equation ([T| gives an estimate of / = 58.9 cm 
using values of uq and a calculated in the effective medium approximation.^ 

A close inspection of Fig. |3](a) reveals that the data above x > 62 cm is actually composed 
of two superimposed peaks that both move to larg Gr X clS db function of d s . Figure Kmc) 
shows two examples of a double-gaussian fit to the data in this region using a standard 
unconstrained, nonlinear optimization routine. The gaussians resulting from the fit are 
shown individually (blue and red) in addition to the combined fit (black). Although the fit 
equation j l e~^ 1 ^ x ~ dp1 ^ 2 + j 2 e~^ 2 ^ x ~ dp2 ^ 2 contains six free parameters, the purpose of the fit is 
to obtain an estimate of the peak positions d p \ t 2 and their relative amplitudes 7^2. Fig. |3](c) 
demonstrates that the data is well described by the double-gaussian, with a low-amplitude 
peak (Peak 1) closer to the GIL and a larger-amplitude peak (Peak 2) farther away. In 
both cases the amplitude of Peak 2 is about three times larger than Peak 1, suggesting that 
Peak 2 is the main focusing peak of the GIL. The relative amplitudes of Peaks 1 and 2 are 
observed to follow the same qualitative behavior for all the source positions in Fig. [3]^a) . 

Figure^d) plots the inverse positions l/d p i t2 extracted from the gaussian fits as a function 
of l/d s . An ideal lens will produce a linear trend with a slope of —1 and an intercept of 
1/f- Although the trends for both peaks are linear and have intercepts that yield similar 
focal lengths, the slope of Peak 1 is much less than that of an ideal lens. However, the 
trend for Peak 2 results in a slope of —1 and its focal length / = 59.3 ± 1.5 cm agrees with 
the estimate of / calculated using Eqn. 0. The dashed line in Fig. |3](d) plots the peak 
locations obtained from the MST calculations in Fig. [3](b) . The theory produces a slope of 
—1 and focal length / = 61.1 cm that closely match both the measured data and the ideal 
lens equations. 



We propose that Peak 1 and the other low-amplitude peaks in Fig. |3^a) are the result 
of constructive interference between waves scattered from the support structure of the lens. 
Low-amplitude, circular interference fringes can be observed in Fig. [2]^a) emanating from 
above and below the plot area centered at x ~ 25 cm. These fringes are the result of scat- 
tering off of stabilizing pillars at the corners of the GIL support structure. While averaging 
over a few initial pulse cycles will reduce the interference, a small number of cycles gives a 
poor approximation to a continuous wave measurement and limits the number of multiple 
scattering events that contribute to the focusing peak. Therefore we have chosen to average 
over many cycles and rely on the gaussian fitting routine to remove the spurious interference. 

Figure [4] demonstrates that our GIL design acts as a lens with the source off the central 
axis. Figure[4|a) shows (P/Po) 2 measured with the spherical source located at a 14.7° angle 
with respect to the origin. Figure W[h) shows the MST calculation for the same source 
location. Thin white lines in Figs. Wa,b) indicate a 14.7° angle with respect to the x-axis 
and extend to the expected focusing positions of an ideal lens with / w 60 cm. Both the 
measured data and the MST calculation demonstrate a strong focusing peak at the expected 
location. Interference fringes from the support pillar can also be observed superimposed on 
the focusing peak in Fig. [4ja). 

In summary, we have designed and constructed a gradient index lens that operates in 
water at sonic frequencies. Our transmission measurements demonstrate that our GIL design 
focuses as an ideal lens based on the optical GIL equations. Our measurements are also 
consistent with the focusing positions obtained from two-dimensional models using multiple 
scattering theory. We emphasize that our GIL behaves as an ideal lens at the limit of 
homogenization (A ~ 4a) and with a thickness on the order of a wavelength (t = 5A/4). 
Such performance at the limit of homogenization theory demonstrates the versatility of 
phononic crystals designed using multiple scattering theory. 

This work was supported by the U.S. Office of Naval Research. 
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FIG. 1. (a) Plan schematic of the gradient index lens, (b) Cylinder radius R plotted vs position 
along the y-axis. (c) Digital photograph of the GIL in the isolation tank. 
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FIG. 2. (a) Normalized pressure amplitude P/Pq plotted vs x and y, measured 2.144 ms after the 
initial pulse leaves the source, (b) Measured, normalized pressure intensity {P/Pq) 2 plotted vs x 
and y after averaging over a 10-cycle pulse, (c) Normalized pressure intensity (P/Pq) 2 calculated 
using MST. The focusing peak maximum is marked by a +. 
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FIG. 3. (a) Measured, normalized pressure intensity vs x for source positions d s = 196, 185.8, 
175.7, 165.5, 155.4, 145.2, 135, 124.9, and 114.7 cm. (b) Normalized pressure intensity calculated 
using MST for source positions d s = 109a, 103a, 98a, 92a, 86a, 81a, and 75a. Inset: expanded 
region of the y-axis showing the focusing peak positions, (c) Two focusing peaks from panel (a) are 
replotted as circles (upper region offset for clarity). Black lines indicate a double-gaussian fit, with 
blue (Peak 1) and red (Peak 2) lines showing the component gaussians individually, (d) Inverse 
peak positions l/d p i,2 plotted vs inverse source positions l/d s . Blue and red lines are fits to the 
trends of Peaks 1 and 2 respectively. The dashed line plots the peak positions of the simulated 
data in panel (b). 
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FIG. 4. (a) Measured, normalized pressure intensity (P/Po) 2 plotted vs x and y with the source 
at a 14.7° angle with respect to the x-axis. (b) Normalized pressure intensity {P/Pq) 2 calculated 
using multiple scattering theory with the source ctt £l 15 ° angle. White lines indicate the off-axis 
angle; positions of the expected focusing peaks are marked with a +. 
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